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Abstract 


Organisms can adapt to habitat conditions that ensure their survival. 
Habitat separation can lead to different populations of body shape during 
the phenotypic plasticity process. Both traditional and geometric 
morphology are being used in fish population studies. In this study, the 
body shape differences between Garra rufa (Heckel) populations 
captured in the Jarrahi River (from the Tigris Basin) were investigated 
using traditional and geometric morphometric methods. The samples 
were captured from the Rostam Abad, Aghajari and Behbahan tributaries 
and transferred to the laboratory. For the traditional morphometric 

Received: 11 March 2020 analysis, 10 meristic characters and 19 morphometric characters were 

pee eo) May eons measured. Geometric morphological information was extracted using 13 

Debra a eg landmark points on left side photographs of individual fish. According to 
the results of the traditional morphometric analysis, there were 
differences between the three populations in meristic (lateral line scales, 
predorsal scales, circamucaudal scales) and morphometric (14 of 19 
characters) traits. In the geometric morphometric analysis, the major part 
of the shape variation is due to landmark points in the head region and the 
dorsal fin base, with the anal fin and caudal peduncle being the most 
conservative body regions. The populations had significant differences in 
body shape with populations from Aghajari and Behbahan tributaries 
being most similar and the Rostam Abad population was different from 
the two other populations. 
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Introduction 


Organisms can adapt to environmental conditions that ensure their survival. During the 
anatomical development history and ontogeny of each organism, the formation of its body 
shape 1s strongly influenced by the environmental conditions, and fish, like other organisms 
can show a similar ability (Nacua et al., 2010). Fish experience different environmental 
conditions throughout their life that can affect various traits including their body shape (Costa 
and Cataudella, 2007). 
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The study of the phenotypic plasticity in populations that live in different environmental 
conditions allows us to better understand changes in population characteristics in relation to 
environmental changes (Kuliev, 1988). Habitat separation can lead to the formation of 
different populations of fish species during the phenotypic plasticity process, due to various 
environmental factors (Wootton, 1990). Morphological characteristics of fish are very 
sensitive to environmental factors such as substrate type, water flow, vegetation, competition, 
predation and availability of food sources (Nicieza, 1995). 


Traditional morphometric approaches like morphometric and meristic characteristics and 
geometric morphometric approaches like landmark-based and truss network have a 
widespread application in fish population investigations, systematics, and separation of fish 
populations (Bakhoum, 1994). These features are used as an effective way to identify 
differentiation and overlap in populations and are the first step in assessing the population 
structure of any species (Turan, 1999). 


The genus Garra Hamilton belongs to the carp family (Cyprinidae) and has at least 100 
species world-wide (Fricke et al., 2019) and 11 species in Iran (Esmaeili et al., 2018; 
Mousavi-Sabet et al., 2019). The species, Garra rufa (Heckel), is a small-sized, tropical and 
widespread freshwater fish that lives in Syria, Iran, Iraq and Turkey (Fricke et al., 2019). 
Previous studies on G. rufa are limited and focused on length-weight and length-length 
relationships (Esmaeili and Ebrahimi, 2006; Gerami et al., 2013; Hamidan and Britton, 2013; 
Hashemzadeh Segherloo et al., 2015; Keivany et al., 2015b; Keivany and Zamani-Faradonbe, 
2017), karyotype (Gorshkova et al., 2012), microbiological risks to humans (Heistinger et al., 
2011), reproductive biology and age determination (Abedi et al., 2011), and mortality (Majtan 
et al., 2012). Some of them also used morphological and morphometric characters to separate 
G. rufa populations (Keivany et al., 2015a; Cicek et al., 2016; Yedier et al., 2016). 


The aim of this study is to investigate the body shape diversity in Garra rufa population(s) in 
the Jarrahi River using two methods, geometric morphometric and traditional morphometric 
methods. 


Material and Methods 


For this study, 55 specimens of Garra rufa were collected using electrofishing techniques 
during October 2017 from three different populations (19 specimens from Rostam Abad, 22 
specimens from Aghajari and 14 specimens from Behbahan tributaries) from the Jarrahi River 
(Fig. 1). The ecological parameters of these sites are in Table 1. Captured fish, after 
anesthesia with a standard solution (1%) of Clove oil, were preserved in 10% neutralized 
formalin solution and transferred to the laboratory for analysis. 


Traditional Morphometric Analysis 


Ten meristic characters including the number of soft rays in the dorsal, caudal, anal, ventral 
and pectoral fins, the number of scales on the lateral line, above the lateral line, below the 
lateral line (to pelvic-fin origin and anal fin origin), circamucaudal and predorsal were 
counted. Also, 19 morphometric characters (Fig. 2) were measured with a digital caliper (0.1 
mm accuracy). The methods for ray and scale counts and measurements follow Kottelat and 
Freyhof (2007). 


The three measurements in the head area, including the snout length, the eye diameter, and the 
postorbital length, were modified by dividing by head length and the other measurements, 
excluding the total, fork and standard length were modified by dividing by standard length. 
Therefore, total length, fork length and standard length were used in the analysis without 
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modification. After modification of the above-mentioned data, all the meristic and 
morphometric characters were analyzed for normality distribution using the Kolmogorov— 
Smirnov test. Nonparametric data were analyzed using the Kruskal-Wallis test and parametric 
data were analyzed using the ANOVA test in SPSS software version 19. 
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Figure 1: Location of the study area in Khuzestan Province and Iran. Map created in ArcGIS 
(version 10.3) with Imagery Basemap. 
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Table 1: Description of some ecological parameters of the sampled sites in this study. 


Population No. Temperature EC TDS pH Habitat Current 
Rostam Abad 19 25.3 3646 1825 8.40 River Moderate 
Aghayjari De WAS) 2582 1291 8.46 Large river Fast 
Behbahan 14 25.8 2484 81242 = 7.93 River Moderate 


No. = Number of fish sampled; EC = Electrical conductivity; TDS = Total dissolved solids. 





Figure 2: The localities of measurements on the left side of Garra rufa from the Jarrahi 
River. 1: Total length, 2: Fork length, 3: Standard length, 4: Predorsal fin length, 5: Dorsal fin 
base length, 6: Post-dorsal fin length, 7: Preanal fin length, 8: Anal fin base length, 9: Post- 
anal fin length, 10: Prepectoral fin length, 11: Prepelvic fin length, 12: Pectoral-pelvic 
distance length, 13: Head length, 14: Snout length, 15: Eye diameter, 16: Postorbital length, 
17: Head depth, 18: Body height, 19: Caudal peduncle depth. 


Geometric Morphometric Analysis 


To extract geometric morphometric information, the left side of the sample fish was 
photographed with a Canon digital camera (model EOS 800D), then 13 landmark points (Fig. 
3) were selected and digitized by the TpsDig2 software version 2.10 (Rohlf, 2010). Then, to 
remove non-shape variations (like location, orientation and scale) from the obtained landmark 
configurations, they were superimposed by GPA (General Procrustes Analysis) (Rohlf and 
Bookstein, 1990) using MorphoJ software version 1.06d (Klingenberg, 2011) and PAST 
software version 2.17b (Hammer et al., 2001). Finally, the Principle Component Analysis 
(PCA), Canonical Variate Analysis (CVA) and Cluster Analysis (CA) were also conducted 
using MorphoJ software version 1.06d and PAST software version 2.17b. 





Figure 3: The localities of the 13 landmarks on left side of G. rufa specimens from the Jarrahi 
River. 
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Results 


Traditional Morphometric Analysis 


Based on the Kolmogorov—Smirnov test, seven meristic characters were nonparametric and 
the results of Kruskal-Wallis analysis of these characters showed that there were no 
differences between the studied populations of G. rufa (Table 2) (p>0.05). Four meristic 
characters were parametric and analysis by ANOVA test showed that the studied populations 
were different in three of them - the number of scales on the lateral line, the number of 
predorsal scales and the number of circamucaudal scales (p<0.05) (Table 3). 


Based on the normality test, all modified and unmodified morphometric character data had 
normal distribution, so they were analyzed with the ANOVA test. The ANOVA test and 
Duncan grouping results are shown in Table 4 and populations were different in many 
characters (p<0.05) except for the preanal fin length, anal fin base length, pectoral-pelvic 
distance length, head depth, and caudal peduncle depth (p>0.05) (Table 4). 


Geometric Morphometric Analysis 


The PCA based on the landmark-based data shows that the first two principal components 
explained 56.21% of shape variation (Fig. 4): PCl = 35.195%, PC2 = 21.015%. The 
populations are therefore partially grouped along the first two axes (PCI and PC2) (Fig. 5). 
Morphological differences between the three populations sampled are visualized through thin- 
plate spline transformation grids and photographic representations using principle axis values 
from the PCA. The major part of the shape variation is due to the landmarks of the head 
region (landmarks 1, 3, 5, 12 and 13) and the dorsal fin base (landmarks 5 and 6), with the 
anal fin and caudal peduncle being the most conservative body regions. Shape deformation 
along PC1 is mostly found on the head region and the beginning and end of the dorsal fin and 
shape deformation along PC2 is mostly found on the head region, the beginning and end of 
the dorsal fin and the beginning and end of the anal fin (Fig. 6). 


Table 2: The Mean + SD and Kruskal—Wallis analysis results for nonparametric meristic characters in 
Garra rufa populations from the Jarrahi River (p>0.05). 
Below the Below the 


; Above the : f Dorsal Anal Caudal Ventral fin 
Populations ene lateral line to lateral line to face tin earn ee 
pectoral fin the anal fin y y y y 
Rostam Abad = 3..334.52 3+0 3+0 8+0 5+0 17+0 7.54.55 
Aghayari oo 3.28+.46 3,.22+.43 8+0 5+0 16.86+.35 7.68+.57 
Behbahan 3+0 3+0 3+0 8+0 5+0 1740 = 7.75+.46 
Pp 0.149 0.155 0.216 1.00 1.00 0.347 0.345 


Table 3: Mean + SD and ANOVA analysis results for parametric meristic characters in Garra 
rufa populations from the Jarrahi River. 


Populations Lateral line scale _ Pectoral"’ Predorsalscales Circamucaudal scales 


Rostam Abad 33.17+40.75° 11.83+41.17 11.5+0.55° 14.67+41.21° 
Aghajari 33.28+0.96° 11.91+40.81 10.8+1.23° 14.56+0.99° 
Behbahan 31.75+1.04" 12.25+0.88 9.12+0.35" 13.7541.19° 


a>. Duncan grouping in ANOVA (p<0.05), "*: not significant (p> 0.05). 
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Figure 4: The percentages of total variance for principal components in Garra rufa 
populations from the Jarrahi River. 
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Figure 5: Scatter plots of Principal Component Analysis (PCA) results for Garra rufa in the 
Jarrahi River. 





Figure 6: Visualized thin-plate spline transformation grids of Garra rufa body shape rather 
than the consensus body shape through PC 1| and PC 2 in the Jarrahi River. 
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In order to assess small population differences in the shape, CVA was performed between 
groups. Based on Figure 7, the populations are strongly grouped along the first two axes (CV 
1 and 2) and these three populations were different in body shape (Wilks lambda= 0.02532 
and p-value = 0.00000005275) (Fig. 7). Morphological differences between the three 
population samples were visualized through thin-plate spline transformation grids and 
photographic representations using canonical axis values from the MANCOVA. 


Shape deformation along CV | is mostly found on the head region and the beginning and end 
of the dorsal fin and along CV 2 it 1s mostly found on the head region, the beginning and end 
of the dorsal fin and the beginning and end of the anal fin. The Behbahan tributary samples 
had bigger heads (change in position of landmarks 1, 3, 4, 12 and 13), shorter dorsal fin base 
lengths (change in position of landmarks 5 and 6), shorter anal fin base lengths (change in 
position of landmarks 9 and 10) and more shallow caudal peduncles (change in position of 
landmark 9) than the Rostam Abad tributary samples (Fig. 8). 
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Figure 7: Scatter plots of Canonical Variance Analysis (CVA) result of the studied 
populations of Garra rufa in the Jarrahi River. 








Figure 8: Visualized thin-plate spline transformation grids of Garra rufa rather than 
consensus body shape through CV 1 and CV 2 in the Jarrahi River. 
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The Mahalanobis distances and Procrustes distances results (Table 5) show that there are 
differences between the body shape of specimens from the Rostam Abad tributary compared 
to Aghajari and Behbahan tributaries, and the differences were significant (p<0.0001). Most 
differences were between specimens from the Behbahan tributary and the Rostam Abad 
tributary. 


The Results in Figure 9 show that the populations of Aghajari and Behbahan tributaries were 
grouped and that the Rostam Abad population was different from the other two populations. 


Table 5: Results of Mahalanobis and Procrustes distances from CVA/MANOVA of Garra rufa 
populations from the Jarrahi River. 








Mahalanobis distances Procrustes distances 
Populations Rostam Abad  Aghajari Populations Rostam Abad _  Aghayjari 
Aghayjari 3.78 Aghayjari 0.043 
Behbahan 4.54 2.81 Behbahan 0.051 0.035 
Distance 
eee é ati 
0.048 0,042 0.036 0.030 0.024 0.018 0.012 0.006 0.0 r-ee ie iii | 
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Figure 9: Distance tree of the Cluster Analyses of the three sampled populations of Garra 
rufa in the Jarrahi River. 


Discussion 


The traditional and geometric (landmark-based) morphometric data of Garra rufa specimens 
used in this study led to the identification of different morphological types. Many fish species 
exhibit morphological differences between habitats (Schluter and McPhail, 1992). The results 
of this study show that there were differences between the studied populations in some 
traditional morphometric traits including meristic characters (number of lateral line scales, 
predorsal scales and circamucaudal scales), morphometric characters (total length, fork 
length, standard length, predorsal fin length, dorsal fin base length, post dorsal fin length, post 
anal fin length, prepectoral fin length, prepelvic fin length, head length, snout length, eye 
diameter, postorbital length, body height) and in results of the geometric morphometric 
analysis. The three populations show differences between them and two populations, the 
Aghajari and Behbahan sites, group together while the Rostam Abad tributary population is 
separate. 


Divergent selection for fish inhabiting different habitats might arise from two major 
differences between these environments: (1) water velocity and (2) resource 
availability/composition (Langerhans et al., 2003). All morphological variations between fish 
populations may be affected by the genetic makeup of individuals, but they may also reflect a 
series of adaptations to a mixture of environmental factors like temperature, salinity, 
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dissolved oxygen, radiation, water depth, current flow influence, food availability, feeding 
mode and habitat use (Antonucci et al., 2009; Fischer-Rosseau et al., 2009; Spoljaric and 
Reimchen, 2011). For example, the Rostam Abad site had more vegetation cover, and more 
electrical conductivity (EC) and total dissolved solids (TDS) than the other sites, but the 
Aghayjari station had more water depth and velocity than other sites. But any discussion about 
the relationship between body shape and habitat characters needs more detailed measurements 
of habitat characters and the appropriate analysis to find these relationships. 


The CVA result shows that the distinction between the Rostam Abad and Aghajari, Rostam 
Abad and Behbahan, and Aghajari and Behbahan tributary populations is comparable using 
Mahalanobis and Procrustes distances showing morphological similarities. In fact, one of the 
valid ways to statistically evaluate the morphological differences between groups 1s the use of 
Mahalanobis and Procrustes distances; although there is a higher Mahalanobis distance 
between the Behbahan and Rostam Abad populations than between the Behbahan and 
Aghajari or Rostam Abad and Aghajari tributary populations. 


The phenotype of an organisms is always interacting with environmental and biological 
factors that directly or indirectly affect the body shape and with selection the favorable traits 
in an environment, increases the chance of survival and reproduction. This can drive the 
formation of a new population and sometimes even speciation (Spoljaric and Reimchen, 
2007). 


Five important morphological measurements in ecology and morphology (Ecomorphology) 
studies of cyprinid species are the position of the dorsal and anal fins relative to the body, the 
position of the mouth, and the position of the dorsal and anal fins relative to one another 
(Wood and Bain, 1995). Fish morphology has been linked to habitat use (Videler and Wardle, 
1991) and body shape can indicate performance of important tasks such as swimming 
(Triantafyllou et al., 2000) and feeding (Nilsson and Broenmark, 2000). Body shape may also 
limit habitat use through biological energy (bioenergetics) constraints. Different habitat types 
become more costly to use because of design limitations (Facey and Grossman, 1990). 
Overall, morphology provides physical and biological energy thresholds that contribute to 
habitat selection and a mechanistic understanding of morphology-habitat associations may be 
used to identify habitat preferences of a species (Schoener, 1986). 


Traditional morphometric and geometric morphometric analysis have many _ benefits; 
including they are inexpensive, all the software used is available as freeware, the techniques 
are faster and easier than molecular techniques, and also are suitable for amateurs (Grinang et 
al., 2019). 


Cicek et al. (2016) studied two species of Garra (G. variabilis (Heckel) and G. rufa) from the 
Tigris Basin, Turkey, using geometric morphometric methods and their results showed that 
there is no difference between these two species in terms of size (centroid size); however, 
they are quite different from each other in terms of shape; and this difference is significant. 
Cicek et al. (2016) suggested that these differences are thought to be caused by feeding habits 
and habitat structure. 


Keivany et al. (2015a) studied 29 populations of G. rufa from six river systems and basins in 
Iran. They measured 28 morphometric characters and counted 10 meristic characters. Among 
morphometric traits, the ratio of the dorsal fin base/standard length and pectoral fin/standard 
length, showed significant differences (p<0.05) and among the meristic traits, except for the 
number of vertebrae, caudal fin rays and pharyngeal teeth, all other traits were significantly 
different (p<0.05). Furthermore, thirteen populations of G. rufa from the Tigris and Persian 
Gulf basins in Iran show significant differences in 28 morphometric and nine meristic 
characters (p<0.05) (Ghalenoei et al., 2010). 
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The effects of environmental conditions are predominant during the early stages of fish life 
and individuals are more sensitive to changes in environmental conditions then (Langerhans 
et al., 2003). Typically, fish that experience the same environmental conditions in their early 
life history have the same body shape (Pinheiro et al., 2005). On the other hand, when the fish 
experience new environmental conditions, it 1s possible for rapid morphological changes to 
occur (Poulet et al., 2004). 


Recent research has shown that morphological differences between different populations of 
the same species can be due to genetic separation and also due to environmental differences 
(temperature, turbidity, food availability, water depth and flow) (Turan, 2000) or 
contributions from all these mechanisms (Brinsmead and Fox, 2002). To determine the effects 
of environmental factors, there 1s a need to study relationships between ecological parameters 
and morphological traits and to determine the contribution of genetics in differentiation of the 
morphology. More studies based on the relationship between molecular characters and body 
shape are also required. 
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